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Fifty-two bats captured during July 2008 in the Philip¬ 
pines were tested by reverse transcription-PCR to detect 
bat coronavirus (CoV) RNA. The overall prevalence of vi¬ 
rus RNA was 55.8%. We found 2 groups of sequences that 
belonged to group 1 (genus Alphacoronavirus) and group 
2 (genus Betacoronavirus) CoVs. Phylogenetic analysis 
of the RNA-dependent RNA polymerase gene showed 
that groups 1 and 2 CoVs were similar to Bat-CoV/China/ 
A515/2005 (95% nt sequence identity) and Bat-CoV/ 
HKU9-1/China/2007 (83% identity), respectively. To propa¬ 
gate group 2 CoVs obtained from a lesser dog-faced fruit bat 
(Cynopterus brachyotis), we administered intestine samples 
orally to Leschenault rousette bats ( Rousettus leschenaulti) 
maintained in our laboratory. After virus replication in the 
bats was confirmed, an additional passage of the virus was 
made in Leschenault rousette bats, and bat pathogenesis 
was investigated. Fruit bats infected with virus did not show 
clinical signs of infection. 

S evere acute respiratory syndrome (SARS) coronavirus 
(SARS-CoV) is a newly emerged zoonotic CoV that 
caused an international epidemic in 2003. Epidemiologic 
studies have demonstrated that the first human cases of 
SARS were caused by CoVs closely related to those found 
in Himalayan palm civets and raccoon dogs in wildlife 
markets (7). This finding accelerated surveys of CoVs spe¬ 
cific for various animals in Southeast Asia to identify reser- 
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voirs for SARS-CoV. These survey findings suggested that 
palm civets and raccoon dogs are an intermediate host of, 
but not a primary reservoir for, SARS-CoV because of the 
low prevalence of SARS-like CoVs in these animals (2). 
Moreover, a large variety of novel CoVs in these surveys, 
including bat SARS-like CoVs, were detected in many bat 
species in the People’s Republic of China and Hong Kong 
Special Administrative Region (5-6). 

Phylogenetic analysis of bat CoVs and other known 
CoVs suggested that the progenitor of SARS-CoV and all 
other CoVs in other animal hosts originated in bats (5,7). 
Recently, bat CoVs in North and South America, Europe, 
and Africa were also reported (8-12). Although extensive 
bat surveys have been conducted, no infectious bat CoVs 
have been isolated from cell cultures, which hinders char¬ 
acterization of bat CoVs and evaluation of the risks posed 
by these viruses to public health. 

In this study, we detected bat CoVs in the Philippines. 
We attempted to isolate bat CoVs and virus RNA from cell 
cultures and from Leschenault rousette bats (Rousettus le¬ 
schenaulti) orally infected with intestinal tissues and con¬ 
tents from a lesser dog-faced fruit bat (Cynopterus brachy¬ 
otis). After infection, clinical signs of infected bats were 
examined, and pathogenesis in bats was investigated. 

Materials and Methods 

Bat Collection 

We obtained 52 bats of 6 species during July 2008 
from Diliman and Los Banos, the Philippines, after receiv¬ 
ing permission from the government. All captured bats 
were anesthetized with an intraperitoneal injection (15 mg/ 
kg) of tiletamine and zolazepam (Virbac, Carros, France) 
and killed by cardiac exsanguination. The experiment was 
conducted in accordance with the Guidelines for the Care 
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and Use of Laboratory Animals, Graduate School of Agri¬ 
culture and Life Sciences, University of Tokyo. 

Extraction of RNA and Reverse Transcription 

Virus RNA was extracted from samples obtained 
from field bats and from experimentally infected bats by 
using an SV Total RNA Isolation System Kit (Promega, 
Madison, WI, USA) according to the manufacturer’s in¬ 
structions. Extracted RNA was eluted in 50 pL of RNase- 
free water. For cDNA synthesis, RNA (5 pL), a random 
hexamer, and a Superscript III Kit (Invitrogen, Carlsbad, 
CA, USA) were used. 

PCR and DNA Sequencing 

All cDNA samples obtained from field bats and exper¬ 
imentally infected bats were tested by using conventional 
and nested PCR. On the basis of previous reports, we used 
a PCR and a pair of consensus primers specific for a highly 
conserved region of the RNA-dependent RNA polymerase 
(RdRp) gene (13). 

Two microliters of cDNA was added to a 25- 
pL reaction mixture containing 2x GoTaq PCR 
Master Mix (Promega) and 0.2 pM of 5'-GGT 
TGGGACTATCCTAAGTGTGA-3' (primer 1) and 5'- 
CC AT CATCAGAT AGA AT CAT CAT A-3' (primer 2). 
The PCR conditions were 2 min at 94°C; 35 cycles for 
20 s at 94°C, 30 s at 50°C, and 30 s at 72°C; and 1 min 
at 72°C. PCR amplicons were gel purified by using Nu- 
cleoSpin Extract II (Machrey-Nagel, Diiren, Germany) 
and cloned by using a TOPO-TA pCR2.1 Cloning Kit 
(Invitrogen). Sequencing was performed in an AB1 3130 
XL DNA analyzer (Applied Biosystems, Foster City, CA, 
USA). 

On the basis of the sequences obtained, we de¬ 
signed new specific primer pairs for the Bat-CoV/Phil- 
ippines/Dilimanl552Gl/2008 sequence (5'-TGATTT 
CT GCAAT GAT ACTT GGTT C-3' and 5ACTT GAT GAT 
CT GT AAC AAC AAT CG-3') and for the Bat-CoV/ 
Philippines/Diliman 1525G2/2008 sequence (5'-TACAAC 
CTACGCTGCAACTC-3' and 5'-ATGAGTGTGCACAA 
GTGCTTAG-3'). These primers were used as the inner 
primer set for the nested PCR after the first PCR was per¬ 
formed with primers 1 and 2. Aliquots (2 pL) of cDNA 
for primary amplification were added to 2x GoTaq Mas¬ 
ter Mix (Promega) and primers 1 and 2. Amplification 
was performed by using 15 cycles at conditions described 
above. Aliquots (2 pL) of primary amplification products 
were used for the second PCR with GoTaq Master Mix and 
the inner primers. The second PCR was performed by using 
35 cycles at the conditions described above. PCR products 
were extracted from gels by using NucleoSpin Extract II 
and subjected to direct sequencing or TA cloning. 


Bat Samples 

Leschenault rousette bats were obtained from zoos in 
Japan. Seven bats were randomly selected for the experi¬ 
ments. In each experiment, 2 bats were placed in a negative- 
pressure isolator. One additional bat was kept in a separate 
isolator as a control. A sample of large intestine from a 
lesser dog-faced fruit bat (C. brachyotis) was homogenized 
in a sterile mortar. After low-speed centrifugation (2,000 x 
g for 10 min), the supernatant was used for oral infection. 
Experimentally infected bats were examined daily for clini¬ 
cal signs of infection. Fecal specimens were obtained from 
a clean translucent plastic sheet spread along the bottom of 
the cage. All bats were killed after being anesthetized with 
diethyl ether, and organs (liver, kidney, spleen, lung, brain, 
and intestine) and serum samples were obtained. 

Detection of Virus mRNA in Bats 

To determine membrane, nucleocapsid, nonstructural 
(Ns)7a, Ns7b, and Ns7c protein nucleotide sequences, we 
conducted PCR and DNA sequencing in the same manner 
as for determination of partial RdRp nucleotide sequence 
described above by using the HKU9-Leader42-64 prim¬ 
er (5 '-CCGTTTCGTCTTGTACGAATCAC-3') and the 
3siteAd20T primer (5 -CT GAT CT AG AGGT ACCGGA 
TCCTTTTTTTTTTTTTTTTTTTT-3'). To detect virus 
mRNA, we conducted reverse transcription-PCR (RT- 
PCR) by using 2 primer sets: HKU9-Leader42-64 and 
N468-448r (5 '-GTT ACGT GT GCCC AT GT CACC-3') and 
HKU9-Leader42-64 and Ns7a440-420r (5'-CAAGCCA 
CAACAACATTAGG-3'). 

Quantitative Real-Time RT-PCR 

cDNA synthesis was performed by using 0.5 pL to¬ 
tal RNA and the PrimeScript RT Reagent Kit (TaKaRa, 
Shiga, Japan) according to the manufacturer’s proto¬ 
col. Virus RNA was quantified by using Power SYBR 
Green PCR Master Mix (Applied Biosystems) with 2 
pL of reverse-transcribed cDNA. Quantitative real-time 
PCR was performed by using the Thermal Cycler Dice 
System (TaKaRa). The temperature program consisted 
of an initial denaturation at 95°C for 10 min, followed 
by 40 cycles at 95°C for 15 s and 60°C for 1 min. The 
primer pair for the real-time PCR was designed on 
the basis of the partial RdRp sequences of Bat-CoV/ 
Philippines/Dilimanl525G2/2008; primers used were 
5 '-T CCT AAGT GT GAT AGAGCT AT GCC-3' and 5'- 
GT GC AC ACT C ATTT GCT AACCG-3'. In each experi¬ 
ment, 10-fold serial dilutions of plasmid DNA containing 
the partial RdRp gene were tested in duplicate to establish 
a standard curve for calculating the relative amount of 
RNA in each sample. All samples were analyzed at least 
3 times. To confirm the specificity of each PCR product, 
we conducted a melting curve analysis immediately after 
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the amplification phase of each PCR. The amount of RNA 
in each sample was expressed as the average value (copy 
number per weight [milligrams] of sample). 

Results 

Virus Detected 

During July 2008, a total of 52 bats were obtained at 
3 locations in Diliman and 1 location in Los Banos, the 
Philippines (Table 1). RT-PCRs for a 440-bp fragment in 
the RdRp gene of Co Vs were performed for large intes¬ 
tine samples, including intestinal contents; 9 (17.3%) of 52 
bats were positive. Differences in the 440-nt sequence in 
the RdRp region were determined after TA cloning of the 
9 positive samples. Sequences indicated that the 2 groups 
of sequences obtained belonged to group 1 CoV (genus 
Alphacoronavirus) (n = 4) and group 2 CoV (genus Beta- 
coronavirus) (n = 5). 

A 440-bp consensus nt sequence of the group 1 
CoV was obtained on the basis of alignment of 4 group 
1 CoV sequences detected (>98% nt identity with each 
other) and deposited in GenBank as Bat-CoV/Philippines/ 
Dilimanl552Gl/2008 (DNA Database of Japan [DDBJ] 
accession no. AB539080). BLAST (www.ncbi.nlm.nih. 
gov/BLAST) search findings of GenBank indicated that 
the partial RdRp sequence was most similar to that of Bat- 
CoV/China/A515/2005 (95% nt identity). 

A 440-bp consensus nt sequence of group 2 CoVs 
was also obtained (>98% nt identity with each other) 
and deposited in GenBank as Bat-CoV/Philippines/ 
Diliman 1525G2/2008 (DDBJ accession no. AB539081). A 
BLAST search suggested that the partial RdRp sequence 
was novel but most similar to that of Bat-CoV/HKU9-l/ 
China/2007 (83% nt identity). A phylogenetic tree was 
constructed with the partial RdRp-deduced amino acid se¬ 
quence (120 aa) and available sequences of known CoVs 
(Figure 1). Data in the tree suggested that Bat-CoV/Phil- 
ippines/Dilimanl552Gl/2008 belonged to group lb CoVs 


and Bat-CoV/Philippines/Dilimanl525G2/2008 belonged 
to group 2d to CoVs. 

Specific and nested primer pairs for group lb bat CoV 
and group 2d bat CoV sequences were designed, and nest¬ 
ed PCR was performed by using cDNAs of all samples. 
Twenty additional amplicons («200-bp sequences) were 
obtained by using primers specific for group 2d bat CoVs. 
After direct sequencing or TA cloning, partial sequences 
of all amplicons obtained were found to be nearly identical 
to group 2d bat CoVs (>98% nt identity) and resulted in a 
total CoV prevalence of 55.8% (Table 1). All sequences of 
group lb bat CoVs were obtained from insectivorous bats 
(4/7, 57.1%), and all sequences of group 2d bat CoVs were 
obtained from frugivorous bats (25/45, 55.6%). 

Virus in Cell Cultures 

Cytopathic effects were not observed in any of the 
cells (Vero E6, Vero, Hrtl8, A549, fcwf-4, BKT-1, Tb-1 
Lu, or primary kidney cells derived from Leschenault rou- 
sette bats) tested with bat intestinal specimens positive for 
both detected viruses by RT-PCR. Results of RT-PCR for 
cell lysates to detect viral replication also were negative. 

Virus Propagation in Fruit Bats 

To obtain bat CoVs from field samples, we admin¬ 
istered virus orally to 2 Leschenault rousette bats main¬ 
tained in the Department of Biomedical Science, Gradu¬ 
ate School of Agricultural and Life Sciences, University 
of Tokyo. The volume of intestine samples collected from 
insectivorous bats was less than that from fruit bats be¬ 
cause of their body size, and all positive samples for the 
group lb bat CoV genome were derived from small in¬ 
sectivorous bats. Oral infection was conducted only with 
samples positive for the group 2d bat CoV genome. A ho¬ 
mogenized large intestine sample (60 mg) derived from 
a lesser dog-faced fruit bat, which contained 7.8 x 10 6 
copies of viral genome, was given orally to 2 fruit bats 
(bats A and B). After confirmation that these bats showed 


Table 1. Prevalence of coronavirus in bats, the Philippines 


Sampling site 

Common name (species) 

No. intestine 
samples tested 

No. positive 
(group 1) 

No. positive 
(group 2) 

Los Banos 

Lesser dog-faced fruit bat (Cynopterus brachyotis) 

4 

0 

2 


Cave nectar bat ( Eonycteris spelaea) 

3 

0 

2 


Greater musky fruit bat ( Ptenochirus jagori) 

14 

0 

11 

Diliman (site A) 

Lesser dog-faced fruit bat ( C. brachyotis) 

1 

0 

1 


Cave nectar bat (E. spelaea) 

1 

0 

1 


Greater musky fruit bat (P. jagori) 

1 

0 

0 

Diliman (site B) 

Cave nectar bat (E. spelaea) 

1 

0 

1 


Java pipistrelle bat ( Pipistrellus javanicas) 

3 

0 

0 


Lesser Asiatic yellow bat ( Scotphilus kuhlii) 

4 

4 

0 

Diliman (site C) 

Lesser dog-faced fruit bat (C. brachyotis) 

18 

0 

6 


Greater musky fruit bat (P. jagori) 

1 

0 

0 


Geoffroy rousette bat ( Rousettus amplexicaudatus) 

1 

0 

1 

Total 


52 

4 

25 
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no clinical signs of infection, they were killed 6 days after 
infection. 

Virus genome was detected only in the small and large 
intestines of both bats by RT-PCR (Table 2). Virus was 
not detected in these intestine samples by cell cultures. 
Virus genome was detected by RT-PCR in fecal samples 
obtained during daily observations for clinical signs, and 


<4 - 

79E 


Bat-CoV HKU6-1 DQ249224 
BatCoV/A620/2005 DQ648828 
Bat-CoV/M.das/Germany/D3 3/2007 EU375854 
Bat-CoV/Philippines/Dilimanl 552G1 /2008 

Bat-CoV/A515/2005 DQ648822 
Baf-CoV/A512/2005 NC 009657 
NC 003436 PEDV 


- Bat-CoV/USA/RM-11 EF544563 

- Bat-CoV/A911/2005 DQ648850 

- BatCoV HKU8-1 DQ249228 

loo r NC 002306 TGEV 
'AY994055FIPV 

-NC 005831 HCOV-NL63 


-EF203064 Bal-CoV HKU2 

I-NC 002645 HCOV-229E 

“Si— Bat-CoV/GhanaKwam/19/2008 FJ710046 
ioo,- NC 001451 IBV 


Group 1 


r£ 


L T urkey co ronavirus AF124991 
so, NC 004718 SARS-human 
AY304488 SARS-civet 
1 NC 009694 Bat SARS CoV HKU3 
— Bat-CoV/GhanaKwam/20/2008 FJ710047 
|— NC 009021 Bat-CoV HKU9 

Bat-CoV/Philiopines/Dilimanl 525G2/2008 

I-NC 009019 Bat-CoV HKU4 

'-NC 009020 Bat-CoV HKU5 

-NC 001846 MHV 
- NC 006577 CoV-HKUl 
NC 005147 HC0V-0C43 
r NC 007732 PHEV 
l NC 003045 BCoV 


Group 3 


Group 2 



Figure 1. Phylogenetic tree based on deduced amino acid sequences 
of partial RNA-dependent RNA polymerase of coronaviruses 
(CoVs), the Philippines. The 2 new viruses detected in this study are 
underlined . Percentage of replicate trees in which the associated 
taxa clustered in the bootstrap test (1,000 replicates) is shown next 
to the branches. The tree is drawn to scale, with branch lengths in 
the same units as those of the evolutionary distances used to infer 
the phylogenetic tree. Evolutionary distances were computed by 
using the Poisson correction method and are shown as number 
of amino acid substitutions per site. All positions containing gaps 
and missing data were eliminated from the dataset. The final 
dataset included 120 positions. Phylogenetic analyses were 
conducted in MEGA4 (14). Coronaviruses used for comparisons 
and their GenBank accession numbers are human coronavirus 
(HCoV) 229E (HCoV-229E) (NC_002645), porcine epidemic 
diarrhea virus (PEDV) (NCJD03436), transmissible gastroenteritis 
virus (TGEV) (NC_002306), feline infectious peritonitis virus 
(FIPV) (AY994055), human coronavirus NL63 (HCoV-NL63) 
(NC_005831), bat-CoV/A512/2005 (NC_009657), bat-CoV/ 

A515/2005 (DQ648822), bat-CoV/A620/2005 (DQ648828), bat- 
CoV/A911/2005 (DQ648850), bat-CoV/GhanaKwan/19/2008 

(FJ710046), bat-CoV/GhanaKwan/20/2008 (FJ710047), bat- 
CoV/M.das/Germany/D3.3/2007 (EU375854), bat-CoV/USA/RM- 
11 (EF544563), bat-CoV HKU2 (EF203064), HKU4 (NC_009019), 
HKU5 (NC_009020), HKU6 (DQ249224), HKU8 (DQ249228), 
HKU9 (NC_009021), CoV-HKUl (NC_006577), human coronavirus 
(HCoV-OC43) (NC_005147), murine hepatitis virus (MHV) 
(NC_001846), bovine coronavirus (BCoV) (NC_003045), porcine 
hemagglutinating encephalomyelitis virus (PHEV) (NC_007732), 
human severe acute respiratory syndrome coronavirus (SARS) 
(SARS-human) (NC_004718), civet SARS-like coronavirus 
(SARS-civet) (AY304488), bat-SARS-like coronavirus HKU3 (bat- 
SARS-CoV HKU3) (NC_009694), infectious bronchitis virus (IBV) 
(NC_001451), and turkey coronavirus (AF124991). 


viral genome copy number was determined by real-time 
RT-PCR (Table 3). Virus copy number peaked on day 3. 
On day 4, a fecal sample was not collected because feces 
were not found on the bottom of the isolator. 

Virus mRNA in Experimentally Infected Bats 

For murine hepatitis virus and several CoVs, an «70- 
bp leader sequence is added to the 5' end of the transcrip¬ 
tion regulatory sequence of each nested mRNA during 
mRNA processing (15,16). For bat CoVFIKU9-l, which 
was most similar to group 2d bat CoV, a complete genome 
sequence and putative transcription regulatory sequence of 
Bat-CoV F1KU9 were predicted (17). On the basis of that 
report, primer HKU9-Leader42-64, including a leader se¬ 
quence, was designed (Figure 2). The FlKU9-Leader42-64 
primer and 3siteAd20T primer, which included the oligo 
dT sequence, were used for PCR with RNA extracted from 
intestines of bats A and B. Amplicons were cloned, and 
partial genomic sequences of group 2d bat CoV mem¬ 
brane, nucleocapsid, Ns7a, Ns7b, and Ns7c genes were 
determined. These sequences were deposited in GenBank 
(DDBJ accession no. AB543561). A phylogenetic tree was 
also constructed with the deduced amino acid sequence 
(463 aa) of the complete N gene of group 2d bat CoV and 
available sequences of known CoVs. The tree showed the 
same topology as that constructed with deduced amino acid 
sequence of the partial RdRP gene. The N gene nucleotide 
sequence was most similar to that of Bat-CoV/FlKU9-l/ 
China/2007 (72% identity). 

To confirm presence of transcribed virus mRNA in 
bats A and B (Table 2), RT-PCR specific for mRNA of the 
group 2d bat CoVs was conducted with FIKU9-Leader42- 
64, N468-448r, and Ns7a44CM20r primers. Virus mRNAs 
were detected in RNA extracted from the small intestines 
of bats A and B (Figure 3). All amplicons were sequenced 
and included the nucleotide sequence of the F1KU9-Lead- 
er42-64 primer sequence at the 5' end of the sequences ob¬ 
tained (Figure 2). These results suggest that virus mRNAs 
were transcribed in bats A and B. 

Experimental Infection of Bats 

To determine whether this 2d bat CoV was pathogenic, 
we experimentally infected 5 R. leschenaulti fruit bats. A 
60-mg sample of the small intestine from bat A was given 
orally to 2 bats (bats C and D), and 500 pL of phosphate- 
buffered saline was given orally to 1 bat (bat E) as a con¬ 
trol. These bats were killed 6 days after infection. Clinical 
signs were not observed in the experimentally infected bats. 
Virus genome amplification was not detected by RT-PCR 
in any samples (serum, brain, kidney, liver, lung, spleen, 
and feces). Flowever, virus RNA was detected in the small 
intestine (Table 4). No pathologic changes were observed 
in the intestines or other organs. 
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Table 2. PCR results for bat coronavirus in fruit bats infected by using bat intestinal samples, the Philippines* 


Bat 

Assay 

Liver 

Kidney 

Lung 

Spleen 

Brain 

Small 

intestine 

Large 

intestine 

Serum 

A 

RT-PCR 

- 

- 

- 

- 

- 

- 

+ 

- 


qRT-PCR 

ND 

ND 

ND 

ND 

ND 

1.25 x 10 6 

3.53 x 10 6 

ND 

B 

RT-PCR 

- 

- 

- 

- 

- 

+ 

+ 

- 


qRT-PCR 

ND 

ND 

ND 

ND 

ND 

1.47 x 10 6 

1.50 xio 6 

- 


*Bat intestinal samples containing 7.8 » 10“ copies of coronavirus genome. Values are copies per milligram. RT-PCR, reverse transcription-PCR; virus 
RNA not detected; +, virus RNA detected; qRT-PCR, quantitative RT-PCR; ND, not done. 


Table 3. Time course of detection of coronavirus viral genome by PCR in feces from 2 fruit bats, the Philippines* 


Days after infection 

Test 0 1 

2 3 

4 

5 

RT-PCR 

- + 

ND 

+ 

Quantitative RT-PCR - - 

5.31 xio 4 1.74 xio 7 

ND 

1.5 x 10 6 

'Values are copies per milligram. RT-PCR, reverse transcription-PCR; 

virus RNA not detected; +, virus RNA detected; ND, not done. 



Because virus growth in fruit bats was weaker than 
virus growth in bats after primary infection with field sam¬ 
ples, experimental conditions were changed. Samples (300 
mg) from the small intestine of 2 bats (A and B) were given 
orally to 2 other bats (F and G), which were killed 3 days 
after infection. Six intestinal samples were obtained from 
each bat to determine site specificity of virus growth. The 
entire intestine (duodenum to the large intestine) was divid¬ 
ed into 6 equal parts (1-6). Virus RNA was detected only 
by RT-nested PCR in the small and large intestines (Table 
4), and no pathologic effect was detected in these bats. Vi¬ 
rus RNA was detected in the lower region (parts 3-4 from 
bat F and parts 3-6 from bat G). Flowever, we could not 
determine the specific site of virus replication. 

Discussion 

After the SARS epidemic in 2003, bats were identified 
as carriers of Co Vs in China. Recently, bat Co Vs have also 
been detected in several other regions, including Germany, 
North and South America, and Africa. In the current study, 
we confirmed the presence of 2 CoVs in bats in the Philip¬ 
pines. Our findings suggest that CoV circulation in bats is 
worldwide. Although only 52 bats were tested, CoV RNA 
was present in 55.8% of large intestine samples from these 
bats. Moreover, all bats tested seemed to be healthy. Thus, 
bats may be persistently infected carriers of Co Vs. These 
data are consistent with results of previous reports of CoV 
detection in bats (3-6). 

RNA of group lb bat CoV was detected in 4 (57.1%) 
of 7 insectivorous bats. All bats positive for group lb bat 
CoV RNA belonged to the same species, the Lesser Asiatic 
yellow bat (Scotophilus kuhlii). However, the partial RdRp 
sequence of the virus was most similar to that of Bat-CoV/ 
People’s Republic of China/A515/2005 (95% nt sequence 
identity), which was also detected in bats of the same species 
in the southern China on Hainan Island (5). The Lesser Asi¬ 
atic yellow bat is distributed widely in eastern Asia in Phil¬ 


ippines, Pakistan, Hainan Island, Taiwan, and Borneo (18). 
High similarities of sequences between group lb bat CoV 
and Bat-CoV/China/A515/2005 suggest that these viruses 
are distributed widely in bats enzootic to eastern Asia. 

Although group 1 bat CoV was detected in 1 species 
of insectivorous bats, group 2d bat CoV was detected in 
4 species of frugivorous bats. Five of the 45 frugivorous 
bats were positive by RT-PCR, and an additional 20 were 
positive by RT-nested PCR (prevalence 55.6%). This find¬ 
ing suggests that replication of group 2d bat CoV in the 
intestine is low. 

The complete N sequence of group 2d bat CoV suggests 
that it is a novel virus and most similar to that of Bat-CoV/ 


TRS 

Leader sequence -- 

BatCoV HKU9-1 (42-76] 

BatCoV HKU9-2 (42-76) 

BatCoV HKU9-3 (41-75) 

BatCoV HKU9-4 (40-74) 

HKU9-Leader42-64 : CCGTTTCGTC TTGTACGAAT CAC 


CCGTTTCGTC TTGTACGAAT 
CCGTTTCGTC TTGTACGAAT 
CCGTTTCGTC TTGTACGAAT 
CCGTTTCGTC TTGTACGAAT 


CACTCTTGA 

ACGAAC 

CACTCTTGA 

ACGAAC 

CACTTTTGA 

ACGAAC 

CACTCTTGA 

ACGAAC 


Read 1 : CCGTTTCGTC TTGTACGAAT CACTCTTGA ACGAAC CAATTATG 


Virus genomic 
sequence 


ATGTTGTATA TATTTCCTTC TATATTTGA 


ACGAAC CAATtATG 
N sequence 


Read 2 : CCGTTTCGTC TTGTACGAAT 
Virus genomic 

sequence : AGGTTTTTGC AGGAATGCTT 


CACTCTTGA ACGAAC CATG 


TTTCCCTTA ACGAAC CUG 
NS7a sequence 


Figure 2. Comparison of mRNA sequences of bat coronavirus 
(BatCoV) with viral genomic sequences. Read 1 was obtained 
by using reverse transcription-PCR and HKU9-Leader42-64 
and N468-448r primers. Read 2 was obtained by using HKU9- 
Leader42-64 and Ns7a440-420r primers. Asterisks indicate 
sequence identity for read and virus genome sequences. 
TRS, transcription regulatory sequence; N, nucleocapsid; NS, 
nonstructural. 
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Figure 3. Bat coronavirus/Philippines/Dilliman1525G2/2008 mRNA 
in experimentally infected fruit bats, the Philippines. Reverse 
transcription-PCR results for small intestines of bats A and B. Lane 
M, 100-bp DNAladder; lane-, nontemplate control. 


HKU9-l/China/2007 (77% aa sequence identity). Woo et 
al. (77) detected Bat-CoV/HKU9-l/China/2007 and classi¬ 
fied the viral nucleotide sequence as that of group 2d CoV. 
Our phylogenetic data (Figure 1) suggest that group 2d 
bat CoV and Bat-CoV/HKU9-l/China/2007 belong to the 
same group. 

We attempted to isolate bat Co Vs from several cell 
lines and primary cultured cells. However, virus replica¬ 
tion was not observed, which is consistent with results of a 
previous report (17). No infectious bat CoV has been iso¬ 
lated from cell culture. In the current study, the amount of 
large intestine obtained per bat was <100 mg. Therefore, 
most samples were inadequate for virus isolation, espe¬ 
cially virus-positive samples for group 1 bat CoVs from 
insectivorous bats. 

To obtain sufficient tissue to isolate virus RNA, we 
attempted to infect fruit bats with bat CoV. Although we 
could not obtain bats of the species from which group 2d 
bat CoV was detected in the field survey, we obtained Le- 
schenault rousette bats from zoos in Japan. In addition, 
Bat-CoV/HKU9-l/China/2007, which was most similar to 
group 2d bat CoVs by phylogenetic analysis, was identified 


in this species in Hong Kong (17). This finding indicates 
that fruit bats can be infected with this virus. No signs of 
clinical disease were observed after oral infection with an 
intestine sample derived from a lesser dog-faced fruit bat. 
However, virus RNA was detected in the small and large 
intestines (Table 2), and these intestinal samples contained 
more genome copies than input copies. Furthermore, virus 
RNA was amplified in fecal samples by real-time PCR, 
and viral mRNAs were detected in bats A and B (Figure 
3). These findings indicate that group 2d bat CoVs can be 
orally transmitted to fruit bats and replicate in them. 

Experimental infection was conducted in fruit bats by 
using tissues from virus-infected bats to determine virus 
pathogenicity. However, infected bats showed no signs 
of a pathologic effect, although low levels of virus RNAs 
were detected in the small and large intestines of these bats. 
These findings suggest that fruit bats can be infected with 
bat CoV without showing any signs of infection. However, 
compared with primary infection by field samples obtained 
from C. brachyotis, the level of viral genome amplifica¬ 
tion was low in experimental infection. This finding may 
have been caused by a difference in viral replication in bats 
of different species. In the field survey, partial nucleotide 
sequences of group 2d bat CoVs, were detected in 4 bat 
species. A high prevalence of virus RNA was observed in 
each bat species (Table 1). These findings suggest that the 
group 2d bat CoVs may infect fruit bats of many species. 
The oral infection study showed that CoV is easily trans¬ 
mitted across species. These results, and the fact that most 
reported bat CoV sequences have been detected in several 
bat species (72,77), imply that interspecies transmission in 
bats may be common. 

Further investigation of bat CoV ecology is needed 
to better understand the risk for infection with this vi¬ 
ms. Knowing this risk could help elucidate emergence of 
SARS. Although we demonstrated in vivo propagation of a 
bat CoV, a bat CoV culture system is needed to obtain ad¬ 
ditional information about this virus. 


Table 4. Results of nested and quantitative RT-PCRs of cDNA from bat samples, the Philippines* 



Bat Liver Kidney Lung Spleen Brain Small intestine 

Large intestine 

Serum 


C - - - - 

D - + (ND) + (ND) + (ND) 

E - - - - 


+ (ND) 

+ (6.57 x I0 4 )t 


+ (ND) 


Intestine section 


1 _2_3_4_5_ 6 

F ------- + (ND) + (ND) 

G - - - - ___+ (ND) + (ND) + (ND) + 

_ (5.91 x IQ 4 ) 

*cDNA was synthesized from bat samples obtained after experimental infection (second passage of the group 2d coronavirus in Leschenault rousette 
bats). Values are copies per milligram. Virus load was quantified by real-time PCR. RT-PCR,reverse transcription-PCR; virus RNA not detected; +, 
virus RNA detected; ND, not detected by real-time PCR in RT-PCR-positive samples. 
tResult of nested PCR. 
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